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Abstract 

It is shown that the /i-adaptive mixed finite element method for the 
discretization of eigenvalue clusters of the Laplace operator produces op¬ 
timal convergence rates in terms of nonlinear approximation classes. The 
results are valid for the typical mixed spaces of Raviart-Thomas or Brezzi- 
Douglas-Marini type with arbitrary fixed polynomial degree in two and 
three space dimensions. 


1 Introduction 

The study of optimal convergence rates for adaptive finite element schemes has 
been carried on by several researchers duriirg the last decades in the case of 
source problems (see, e.g., [22, 41, 17, 4, 16, 34]) and more recently has been 
applied to eigenvalue problems as well (see, e.g., [28, 32, 13] for convergence 
aird [20, 14, 19, 12] for optimal rates). Some survey papers are available; we 
refer, in particular, for further reading and references, to [37, 38, 11]. In the 
case of eigenvalue approximation, it has been recently observed that adaptive 
schemes driven by the error indicator associated to an individual eigenvalue may 
produce unsatisfactory results, and that eigenvalues belonging to clusters have 
to be considered simultaneously (see, in particular, [25, 26, 27]). 

In this paper, we study the adaptive approximation of the Laplace eigenvalue 
problem by mixed finite elements. The analysis of the underlying formulation, 
which fits the framework of (0,g)-type mixed problems, is not a mere gener¬ 
alization of the case of standard conforming Galerkin approach (see [6], where 
the convergence and the a priori estimates are recalled). This causes additional 
technical difficulties which were in previous works [24] circumvented by showing 
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equivalence with some nonconforming but elliptic finite element formulation. 
Typically, residual-based a posteriori error estimates are derived by exploiting 
the fact that the error of the eigenvalues as well as the error of the eigenfunctions 
in some weaker norm (usually the norm) is of higher-order compared with 
the error in the energy-like norm. The higher-order convergence, however, is 
not valid in its original format in mixed FEMs, and one technical tool we make 
use of is a fairly abstract superconvergence result for eigenvalue problems where 
a certain error quantity is shown to be of higher order in the norm. For the 
low-order case a similar result was shown in [29] by using the representation in 
terms of nonconforming finite elements from [24]. 

We follow the argument of [17] in order to show the optimality of an adap¬ 
tive finite element scheme which is constructed taking into account clusters of 
eigenvalues in the spirit of [25]. In order to obtain the result, we need to derive 
estimates which are essentially different from the case of standard FEMs: this 
is one of the main contributions of our paper. 

Previous a posteriori estimates for mixed formulation (source or eigenvalues 
problem) mostly showed efficiency and reliability with respect to the vector 
variable only (see [1] and [18, 34]; other results in this context can be found 
in [9, 44, 30, 36, 35]). Estimates involving the scalar variable were present in [24] 
(where, as already mentioned, the equivalence with nonconforming schemes is 
exploited) and in [10] (where the source problem is considered). Another main 
contribution of our analysis is that we show optimality also with respect to 
the scalar variable (see Definitions 6 and 7). This is performed by a suitable 
definition of the error indicator (see Definitions 5 and 9); this allows to prove the 
optimal convergence rate not only for the eigenfunction but for the eigenvalues 
as well (see Section 5). 

The outline of the paper is as follows: Section 2 introduces the problem 
we are dealing with, Section 3 describes the error indicators and our adaptive 
algorithm, Section 4 states the main theorem of our paper, concerning the con¬ 
vergence of the adaptive scheme in terms of a theoretical error indicator which 
is equivalent to the error indicator used for the design of the AFEM algorithm. 
Section 5 shows that the convergence of the error indicator, which is related to 
the convergence of the eigenfunctions, actually implies the convergence of the 
eigenvalues as well. Finally, Section 6 contains all technical results which are 
used in the proof of our main theorem and Section 7 discusses the extension to 
three space dimensions. 

Throughout this paper, we use standard notation for Lebesgue and Sobolev 
spaces and their norms. The norm of a function v over some domain oj is 
denoted by Ijr’ljtj and, if there is no risk of confusion, we write \H = Iklln for 
the physical domain D. The scalar product of LF'{VI) is denoted by (•,•)■ If 
is a disjoint union of subdomains of D, typically a (subset of a) triangulation, 
then 11^1134 = ll'^^llw- We denote the scalar curl of some two-dimensional 

vector field "0 by curl^ = d 2 'ipi — diip 2 and the vector curl of a scalar-valued 
function v by curlw = {—d 2 V,div)"^. In three dimensions we define as usual 
curl'f/j = V X xp. 

The notation A < B refers to an inequality A < CB up to a constant C 
that is independent of the mesh size. We do not trace the explicit dependence 
of the constants on the eigenvalues, cf. Remark 1. 

The mesh-size is typically denoted by /i; when a triangulation Th is obtained 
as a refinement of a given mesh, we denote by Th the coarser mesh. When 
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dealing with the adaptive scheme, we denote by £ the level of refinement, so 
that 7^+1 is the next triangulation in the algorithm obtained from 7^. 


2 Setting of the problem 

Our main result is valid both in two and three dimensions. From now on, 
we discuss the two dimensional setting. Section 7 extends the result in three 
dimensions. 

Given a polygonal domain O, in this paper we are interested in the following 
eigenvalue problem associated with the Laplace operator in mixed form: find 
A G M and u G L?{Vl) with ||u|| = 1 such that for some cr G 77(div; O) it holds 



2.1 Abstract mixed eigenvalue problem 

We cast this problem within the standard setting of abstract eigenvalue problems 
in mixed form of the second type (see [8, 6]). 

Let E, M, Ti, be Hilbert spaces such that M (ZT-L Q M* and consider two 
bilinear and continuous forms a : E x E —?► R symmetric, and & : E x M —>■ R 
which satisfy the usual hypotheses for mixed problems [7]: a is elliptic in the 
kernel of b and b fulfills the inf-sup condition. Moreover, the form a is supposed 
to be positive definite so that the associated norm | • |a is well defined. In the 
pivot space % we consider the scalar product (•,•)■« £md corresponding norm 


In this framework, the continuous eigenvalue problem reads: find A G R and 
u £ M with WuW-H = 1 such that for some cr G E it holds 


a(cr, r) -I- 6(t, m) = 0 Vr G E 
b{u,v) =Vn G M 


( 2 . 1 ) 


and, given finite dimensional subspaces Efi C E and Mh C M (typically associ¬ 
ated to a finite element mesh Th), its discrete counterpart is: find A?i G R and 
Ufi G Mh with ||m;i||h = 1 such that for some ah G E/j it holds 


a{ah,T) + b{T,Uh)=0 'ir £T.h 
b{ah, v) = -Xhiuh, v)n G Mh- 


( 2 . 2 ) 


The following three assumptions ensure the good approximation of the eigen- 
modes (see [8, 6]), where p{h) tends to zero as h goes to zero and Eq and Mq are 
the subspaces of E and M, respectively, containing all solutions to the source 
problem associated with (2.1) when the datum is in 77; the discrete kernel as¬ 
sociated to the bilinear form b is as usual defined as 


K/i = {t G E/i : 6(t, ?;) = 0 Vu G Mh}. 
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Fortid condition. There exists a Fortin operator : Eg —>■ such that 

b{a — nf cr, u) = 0 Vu G Mh 


and 

\(T-'n^FMCr\a< p{h)\\a\\Y.a VctGEo. 

Weak approximability of Mg. 

K'^hjV) < p{h)\Th\a\\v\\Mo Vu G Mg Vt/, G K/,. 

Strong approximability of Mg. 

inf ||u - < p(h)||u||Mo Vu G Mg. 

vh&Mh 

We consider a problem associated with a compact operator, so that the 
eigenvalues are enumerated as 

0 < Ai < A 2 ^ Ag < ... 

(we repeat the eigenvalues according to their multiplicities); the corresponding 
eigenfunctions are denoted by {(cri, ui), (erg, W 2 ), • ■ •} and the {ui}’s form an 
orthonormal system in H. In particular, we have |cri|„ = Xi and HuiU-H = 1 for 

i = 1,2, _We denote by E{X) the span of the {ui}’s corresponding to A. 

Analogously, the discrete eigenvalues can be enumerated as follows 


0 < A/j^ < Xh,2 < • • • < A/i,7v(/t) 

with corresponding eigenfunctions {(cr/i,i,, (cr?j, 7 V(/i),where 
N{h) = dim{Mh) and the {u/t,i}’s form an orthonormal system in H. Here 
we have \ah,i\l = Xh,i and \\uh,i\\n = 1 for z = 1,2,.. .,N{h). 

For a cluster of eigenvalues A„+i,..., A„+i\i of length N G N, we define the 
index set J = {n + 1,..., n + N} and the spaces 

W = span{uj I j G J} and W-r^ = Wh = span{it?i_j | j G J}. 

2.2 Some useful operators 

Definition 1. For any w G M we define G(w) G E as the solution to 

a(G{w),T) + b{T,w) = 0 forallTGT,. (2-3) 

For any Wh G Mh we define its discrete counterpart Gh{wh) G E/j via 

a{Gh{wh),Th)+ b{Th,Wh) =0 for all Th € T,h. (2.4) 

We explicitly notice that when two meshes Th and Th are present, it is important 
to distinguish between Gh and Gh- 

In many applications and corresponding instances of a and b, the above 
definition is related to an integration by parts formula where G(w) is some 
derivative of w. For instance, in the case of mixed Laplacian, G{w) is the 
gradient of w. 
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Definition 2. The solution operators T : TL ^ M and A : ^ S are defined 

by 

{a{Ag,T) + b[T,Tg)=Q Vr € S 

\ b{^9, v) = -{g, v)-H Wv G M 

and Th '.TL ^ Mh and Ah : "H —>■ S/t are their discrete counterparts 

f a{Ahg, Th) + b{Th, Tug) = 0 e ^ 

\h{Ahg,Vh) =-{g,Vh)n Vvh G Mh- 


Definition 3. The operator T^ : TL ^ Mh (^A € Rj is defined by 
( a{Gh{T^g),Th) +h{Th,T^g) =0 Vr/, e 

1 b{Gh{T^g), Vh) = -{Xg, Vh)-H '^Vh G Mh, 


that is, T^ = XTh- 

Let denote the 7t-orthogonal projection onto Wh ■ The following defini¬ 
tion is crucial for the definition of our theoretical error indicator. 


Definition 4. 


The operator Ah : E{X) —> Wh is defined as follows: 


Ah = P}^oT^. 


For the sake of simplicity, we do not include the dependence from A in the 
notation for : it will be clear from the context that when Ah is applied to an 
element of E{X), the corresponding value of A should be used for its definition. 

Lemma 2.1. The operators P^ and T^ commute, that is Ah = P^ = 
T^ o P^. In other words, if {a, u) is an eigenfunction associated with A, then 
AhU solves 


j a{Gh{Ahu),Th) + b{Th, Ahu) =0 Vr/, e Yh 
\ b{Gh{Ahu),Vh) = -{XP^u,Vh)'H ^Vh G Mh. 

Proof. We adapt the result of [27, Lemma 2.2]. The expansion of AhU reads as 
AhU = 'Ylij^j{'Eh'^,Uh,j)'HUh,j, thus AhU solves the discrete linear system (2.6) 
with right-hand side g = J2jej('^h''^^'^h,j)'HXh,jUh,j. For any j G J we have 

Xh,j{ThU,Uh,j)u = -K^h,j,Thu) = a{Gh{Thu),ah,j) = -b{Gh{Thu),Uhj) 

= X{u,uh,j)n, 

which gives the final result that AhU solves the discrete linear system (2.6) with 
right-hand side g = Y.j(ijX{u,Uh,j)HUh,j = XP^u. 

□ 


3 AFEM algorithm and error quantities 

As already mentioned, we are interested in the Laplace eigenvalue problem in 
mixed form with Dirichlet boundary conditions. Namely, with the notation 
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introduced in Section 2, we are making the following choices: 


S = H{div;n) 

M =n = L^{n) 

a{a,T) = (cr,r) 

6(r, v) = (div r, v) 

for an open, bounded, simply-connected polygonal Lipschitz domain fi. 

It follows, in particular that the seminorm | • |a is the norm in Our 

airalysis applies to more general operators (for instance, Neumann boundary 
conditions or non-coirstant coefficients), but we stick to this simpler example 
for the sake of readability. 

We discretize the problem with standard mixed finite elements (including 
Raviart-Thoirras, Brezzi-Douglas-Marini, etc.), see [7] for irrore detail. It is 
well-kirown that this choice satisfies the assumptions discussed in Section 2 
(see, for instance, [8]). 

Moreover, we observe that the following relation (part of the commuting 
diagram) holds true: 

div(S,,) = A/4 (3.1) 

Let us first introduce our error indicator. 

Definition 5. Let Th be a triangulation of LI and let {(Jh,j,Uh,j) G '^h x AL/^ 
he a discrete eigensolution computed on the mesh Th- Then, for all T € Th we 
define 

Vh,j{Tf = WhricThj - Vu/tj)|l| -b \\hT cmlahjWr + hE\\[(Th,j]E ■ 4|||, 

Eee{T) 

where hr is the diameter of T, £{T) denotes the set of edges of T, Le is the 
length of the edge E, and tE is its unit tangent vector. As usual, {<Jh\E • tE 
denotes the jump of the trace of ah ■ tE for internal edges and the trace for 
boundary edges. 

Given a set M. of elements ofTh, we define 

VhAMf = Y. VkATf- 

TeM 


3.1 Adaptive algorithm 

The adaptive algorithm consists of the standard four steps: solve, estimate, 
mark, and refine. In the description of the fours steps, we describe how the 
algorithms runs from level (. to i 1. 

Solve. Given a mesh T the algorithm computes the eigensolutions of (2.2) 
belonging to the cluster {Xej,aej,Uij) for j S J. We assume that the 
discrete solution is computed exactly. 

Estimate. The algorithm computes the local contributions of the error estima¬ 
tor for the eigenfunctions in the cluster {j G J). 
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Figure 3.1: Possible refinements of a triangle T in one level in 2D. The thick 
lines indicate the refinement edges of the sub-triangles as in [5, 42]. 


Mark. The algorithm uses the well known Dorfler marking strategy [22] . Given 
a bulk parameter 9 G (0,1], a minimal subset A4e C Te is identified such 
that 

j&J j&J 

The elements belonging to M.f, are marked for refinement. 

Refine. A new triangulation Tn+i is generated, as the smallest admissible re¬ 
finement of Tt satisfying M-t fl Ti+i = 0 by using the refinement rules 
of [5, 42]. Figure 3.1 shows possible refinements of a triangle. 

To summarize, the adaptive algorithm accepts as input the bulk parameter 
9 and the initial mesh To (with proper initialization of refinement edges as 
in [5, 42]), and returns as output a sequence of meshes {Te} and of discrete 
eigenpairs 

Finally, we shall make use of the following notation: given an initial mesh 
To, regular in the sense of Ciarlet, we denote by T the set of admissible meshes 
in the sense that a mesh in T is a refinement of To obtained using the rules 
of [5, 42]. 

3.2 Error quantities and theoretical error indicator 

The following definition introduces a metric in M. 

Defiuitiou 6. d : M x M is defined as 

d{v,w) = ^\\v - wl|2 -h \G{v) - G(u;)12 

When V (resp. w) belongs to Mh, then Gh{v) (resp. Gh{w)) should be used. 

Remark 1. We note that it may be useful to balance the terms in the square 
root of Definition 6 in terms of X. In particular, if v and w are related to 
eigenfunctions with frequency X, the right scaling would involve multiplying by 
X the term jju — wjj. This is of particular interest it one aims to quantify the 
conditions on the initial mesh-size. In this paper, we do not aim at such a 
quantification and refer the interested reader to [27] for such a X explicit analysis 
in the context of conforming standard finite elements. 

This distance allows us to evaluate the gap between discrete and continuous 
eigenfunctions in the cluster. 

Defiuitiou 7. The following quantity measures how combinations of eigenfunc¬ 
tions in the cluster W are approximated by their discrete counterparts in Wh ■ 

5{W,Wh) = sup inf d{u,Vh). 

ueW -vh&Wh 
\\u\\ = l 
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Given a refinement 7^ G T of the initial mesh To, our theory is based on the 
introduction of the following non-computable error indicator pi which will be 
proved equivalent to the computable indicator rji. 

Definition 8. Let 7^ G T &e a triangulation and for all T £ Th and gu S Mh 
let us consider the following seminorm 

\9h\l^T = WhriGhigh) - '^gh)\\T + ||/iTCurlG/,(5;0||T 

+ ^ hE\\[Gh{gh)]E -tEWl;, 

EeeiT) 

so that 

Vh,j{T) = \uh,j\r^,T- 

Then, given an eigenfunction (a, u) associated to the eigenvalue \ (in particular, 
this is used in the definition of ), we define 

Ph{u-,T) = 

Given a set A4 of elements ofTh, we define 

Ph{u-,Mf = ^ ph{u\T)‘^. 
tgM 

The next lemma is of technical nature and gives a criterion for linear inde¬ 
pendence. It generalizes [13, Prop. 3.2]. 

Lemma 3.1. Recall the notation N = card(7) and suppose that 

e = max \\uj — A/iUj|| < -\/l -|- 1/(2N) — 1. (3.2) 

jeJ 

Then, {AhUj}j^j forms a basis ofWh- For any Wh € Wh with ||w/i|| = 1, the 
coefficients of the representation Wh = are controlled as 

^l7,f<2 + 4N. (3.3) 

jeJ 

Proof. The proof employs Gershgorin’s theorem. Since the proof follows verba¬ 
tim the lines of [27, Lemma 5.1], it is omitted here. □ 

The following lemma states the equivalence between the two introduced esti¬ 
mators. It is clear that the adaptive algorithm will make use of the computable 
indicator g, while the indicator p will be used for the analysis. 

Lemma 3.2 (Local comparison of the error estimators). Provided the initial 
mesh-size is small enough such that (3.2) is satisfied, it holds for any T G Th 
that 

N”i^/r;,(u,;r)2< ("D {2N + LhiupTf 

j&J ^ ' jeJ ^ ^ j&J 


where [A, B] denotes a real interval containing the (continuous and discrete) 
eigenvalue cluster and N is the number of eigenvalues in the cluster. 



Proof. The proof follows from a perturbation analysis as in [27, Prop. 5.1]. 
We include the proof for self-contained reading. Let k G J and consider the 
expansion of AhUk = with coefficients = {AhUk,Uhj). The 

definition of Ah and the symmetry yield 

Ij ~ {AhUkjUhj) = {ThUkjUhj) = iThUk) 

~ ,^h{ThUk)) = Ufc), rt/ij) = Xhj^k{uk,Uh,j)- 

Since {uhj}jeJ is an orthonormal system, we arrive at l] < (B/A)^ 

which implies 

^jeJ ' jeJ ^ ^ jeJ 

This proves the first stated inequality. 

Lemma 3.1 shows that there exist real coefficients {6j | j G J} such that 

Uh.k = E SjAhUj and < 2-h4N. 

leJ jeJ 

The triangle and Cauchy inequalities lead to 

Wh,k\'^^T — f E/ E/ ^ (2 + 4N) - 2 ". 

^jeJ ' 3&J jeJ 

This shows the second stated inequality and concludes the proof. □ 

4 Optimal convergence of the adaptive scheme 

In this section we state the main theorem showing the optimal convergence of 
our adaptive scheme and sketch the principal lines of its proof. The structure 
of the proof is closely related to [17] and relies on several intermediate results 
which, for the sake of readability, will be postponed to Section 6. 

As usual in this context, the convergence is measured by introducing a suit¬ 
able nonlinear approximation class in the spirit of [5]. For any m G N, we denote 

by 

T(m) = {T G T j card(T) — card(7o) < m} 

the set of admissible triangulations in T whose cardinality differs from that of 
To by m or less. 

The best algebraic convergence rate s G (0, -too) obtained by any admissible 
mesh in T is characterized in terms of the following seminorm 

|kFU« = sup to'* inf d{W,W'r). 

meN TeT(m) 

In particular, we have [VF]^,, < oo if the rate of convergence 5{W,W-r) = 
0{m~^) holds true for the optimal triangulations T in T(to). 

The main results of this section, stated in Theorem 4.1, shows that the 
same optimal rate of convergence is reached by the error quantity 5(W,Wj-^) 
associated with the mesh sequence {7)} obtained from the adaptive algorithm 
presented in Section 3. 
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Theorem 4.1. Provided the initial mesh-size and the hulk parameter 9 are small 
enough, if for the eigenvalue cluster W it holds < oo, then the sequence 

of discrete clusters We computed on the mesh Te satisfies the optimal estimate 

S{W,Wi)icaTd{Ti) - card(ro))* < 


Proof. We follow the lines of the proof of Theorem 3.1 in [27]. The main argu¬ 
ments are the same as in [17]. 

Given a positive (3, we consider the quantity 

+/3^d{uj,AeUjf 
jeJ jeJ 

which will be used in the contraction argument of Proposition 6.11. We do not 
consider the trivial case = 0. Choose 0 < r < \W\^_^J^q, and set e{tj = \fT f,i. 
Let N[t) G N be minimal with the property 

It can be easily seen that N{1) > 1, otherwise 

< eW 

but this, together with the definition of would violate the contraction 

property of Proposition 6.11. 

From the minimality of N{P) it turns out that 

N{t) < 2|W|y"e(£)-i/" for all I G Nq. (4.1) 

Let 7^ G T denote the optimal triangulation of cardinality 

card(7^) < card(7o) -I- 

in the sense that the operator A = A.^ of Definition 4 with respect to the mesh 
Te satisfies 

Y^d{u„lu,f < < e{if. (4.2) 

jeJ 

Let us consider the overlay Te, that is the smallest common refinement of Te 
and Te, which is known [17] to satisfy 

card(7£ \Te) < card(7^) — card(7^) < card(7^) — card(7o) < N{i)- (4.3) 

This relation and (4.1)-(4.3) lead to 

card(r^ \ fe) < N{£) < (4.4) 

Let A denote the operator A.^^ with respect to the mesh Te- 
The following estimate 

J2diuj,Aujf < (4.5) 

j&J 
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follows from the quasi-orthogonality (see Proposition 6.9) applied to Th = % 
and Th = Ti- Indeed 

(1 (i7qop(/io ) ) d(Uj , AUj ) ^ (1 ~t“ C^qoP(ho)) E d{uj,Auj)'^. 

j&J jeJ 

Estimate (4.5) follows from the mesh-size condition Cqop{ho) < 1/2 and (4.2). 
We now show the existence of a constant Ci such that 

(4.6) 

j&J j&J 

From the triangle inequality and the discrete reliability (see Proposition 6.7) 
we obtain for any j & J 

d{uj,AiUj)^ < 2d{uj,AiUj)^ -I- 2d{AiUj, AiUj)^ 

< 2d{uj, AiUj)^ -\- 2C^^gj/r^(7/ \ Ti)^ 

+ Cp{ho)^{d{uj, A^Uj) + d{uj,AiUj)Y. 

Provided the initial mesh-size is sufficiently small, this leads to some constant 
C 2 such that with (4.5) it follows 

Y.d{u,,A,u,f <C2e[^f + C2 C'drei pe{uj,Te \ Tg)^. 

j&J 

Let Ceq denote the constant of C 2 Q < Cf.qY^j^jPgiuj.Tg)'^ (which exists 
by reliability). The efficiency (6.2), the definition of e{tj, and the preceding 
estimates prove 

Ceff" E ^ pg{u,,Tg \ fgf 

jeJ jeJ 

< rCeq ^ pg{u,,Tgf + C2CI,, Y. \ 

i6J ieJ 

Defining Ci = {C~g — TCeq)~^C 2 C^^^i, which is positive for a sufficiently small 
choice of r, we obtain (4.6). 

In order to conclude the proof, we now make the following choice for the 
parameter 9: 

0 < 0 < 1 /(C'i(B/A) 2 ( 2 N 2 - 1 - 4N3)) . 

The marking step in the adaptive algorithm selects A4g C 7/ with minimal 
cardinality such that 


jeJ j&J 

Estimate (4.6) and the definition of 9 imply together with Lemma 3.2 that also 
7i\Tg satisfies the bulk criterion, that is 

jeJ jGJ 


II 


The minimality of and (4.4) show that 


card(7W£) < card(7^ \fi) < 


It is proved in [5, 42] that there exists a constant Cedv such that 



<2C'BDv|W|y;T-l/(2«)^4”Vs^ 


The contraction property from Proposition 6.11 implies Q < for k 

0,..., £. Since p2 < 1, a geometric series argument leads to 


i-i 


t-i 





The combination of the above estimates results in 


card(7£) — card(7o) 


<2C'BDv|VP|y;T-fo(2«)^-lApl/(2^) 



The equivalence of with the error d{uj, AiUj)"^ (reliability and efficiency, 

see Section 6) concludes the proof. 


□ 


5 Convergence of eigenvalues 

The previous analysis shows that the adaptive procedure leads to the conver¬ 
gence of the quantity 5{W, Wi) which is related to the eigenfunctions belonging 
to the cluster. In this section we show that this estimate actually implies the 
optimal convergence of the eigenvalues. 

The next discussion has been inspired by [21]. However, we do not make 
use explicitly of the spectral projections and follow a somehow more natural 
argument (at least for symmetric problems). 

As usual, we consider the eigenvalues = If \i (f = 1,...) of T and = 

l/A^^i (i = 1,..., dim(M£)) of and discuss the convergence of to for 

j G J. This standard notation conflicts with our theoretical error indicator; 
nevertheless, we believe that this overlap is not a source of confusion, since it is 
limited to this section where the error indicator is not mentioned. 

Let E : y. y. denote the y projection onto W and E( : y ^ y the y 
projection onto W^. We denote by Fi the restriction of Ei to W 


Fe = Ei\w- 


The following proposition shows that for £ large enough the operator Fi is a 
bijection from W to (which have the same dimension N). 
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Proposition 5.1. For I large enough the operator F( is injective. Moreover, 
Fj~^ is uniformly bounded in C{Wf,,W) and 

sup \\F^^x - x\\n < C5{W, We). 

xeWf 

Proof. It is enough to show that for £ sufficiently large \\Fey — y\\fe < (1/2)||2/||^ 
for aW y & W (see also [21, Lemma 2]). Indeed, from the definition of Fe it is 
immediate to get 

\\Fiy - yWn < \\y - yeWn '^Vi&We 

which implies 

\\Fiy - y\\n < S{W, We)\\y\\n- 

We can then conclude our proof from Theorem 4.1 observing that 5{W,Wi) 
tends to zero. 

□ 

Let us define the following operators from W into itself: 

f = T\w, Te = Ff^TiFi. 

It is clear that the eigenvalues of T {Te, resp.) are equal to pj {pe,j resp.), j G J. 
Lemma 5.2. The following estimates hold true for all x G W 

UT-Te)x\\n<C6{W,We), 
\{A-Ae)x\a<C6{W,We), (5.1) 

\\{A-Ae)x\\s<CS{W,We). 

Proof. Let us denote u = Tx, ue = Tex, a = G(u) = Ax, and ae = Ge{ue) = 
Agx. 

In order to prove the first estimate, we use a standard duality argument and 
introduce the following auxiliary problem: find C G ^ s-nd w G M such that 

J a(C, r) + 6 (t, w ) = 0 Vr G E 

{b{C,v) =-{u - ue,v)-H 'ivGM. 

We clearly have ||CI|i: + I|w^I|m<C'||'« — By standard arguments we get 

||u - = {u-ue,u- ue)n = -KC, u - ue) 

= -b{( - u) - bfllp^eC, u - ue) (5.2) 

= a(cr, C - n_F,fC) + a{G{u) - Ge{ue), n_F,fC)- 

For all ve G Me, the first term can be estimated as follows: 

\a{a, C - LIf.^C)! = l“(o' “ ^e{ve),C - LIf/C) + a{Ge{ve),C - LIf.^C)! 

= |a(cr - Ge{ve),C- IIf.K) - HC - AF,eC,Vi)\ 

= \a{a - Geive),C - Bf/OI 

< Clcr - G£(u£)|a||Clls < C'|cr - Gi{vi)\a\\u - ue\\n- 
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The second term in the last line of (5.2) can be estimated as follows: 
|a(G(w) - n^^X)! < C^|G(u) - G,(u,)U||CIIe <C\{A- Ai)x\a\\u - u,\\n. 

Hence 


\\Tx - TtxWu <C {\(T - Gi{ve)\a + \{A - Af)x\a) Vve £ M^. 

Since the first term is bounded by 5{W^ Wt), the final estimate will follow from 
the second estimate in (5.1). 

Let us prove the second estimate in (5.1). 

From the definition of W we have 

X = OLjUj, 


where we recall that {Xj,aj,Uj) is the generic eigensolution belonging to the 
cluster W and the coefficients are given by aj = {x,Uj). 

Hence, Ax = G(m) with u = Tx and 


jeJ 


1 


^3 


Analogously, from (2.7), 

3^J ^ 


We then obtain 


\Ax - A(,x\a = 





a 


We now show that \aj — G^{T^^Uj)\a can be bounded by S{W,Wg). For all 
£ Ml we have 

\aj - Gi{T^Mj)\l = a{aj - Gi{T^^Uj),aj - Gi{T^Mj)) 

= a{aj - Gi{T^Mj), aj - Gi{ve)) 

+ a{aj - Gi{T^^Uj), Gi{vi) - Gi{T^^Uj)). 

Since the last term is vanishing for the properties of aj and the definitions of 
and G^, we obtain 


kj - Gi{T^^Uj)\a < inf IcTj 


Gi(vi)\a < CSiW,Wi). 


From [7, Prop. 4.3.4] and the definitions of A and Ai it follows that 

||(A - A^)x||e < C|(A - Ai)x\a + C\\x - xi\\h, 

where xi £ Mi is the % projection of x. The first term is readily bounded 
by 6(W, Wi), while the second one is smaller than ||a; — Xa;|X which has been 
already estimated in the proof of Proposition 5.1. 

□ 
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The following proposition is a crucial step for the bound of the eigenvalues. 

Proposition 5.3. The following estimate holds true 

\\f-fi\\ciw)<CS{W,Wif (5.3) 

Proof. Let us define Se = — I : "H —>■ "H. From the boundedness of the 

involved operators, it is immediate to observe that Si is uniformly bounded. 
For all X € W we have 

{f - fi)x = {T - Ti)x + Si{T - Ti)x (5.4) 

since EiSt = 0. Let us estimate the first term. For all x,y G W with ||a:||-H = 
\\y\\n = l 

{{T - Ti)x, y)-H = -h{Ay, (T - Ti)x) + a{{A - Ai)x, Aty) + b{Aiy, (T - Ti)x) 
= -b{{A - Ai)y, (T - Ti)x) + a{{A - Ai)x, Agy). 

The hrst term is bounded by a constant times 5{W^ WiY, while the second 
term can be estimated as follows. 

a((A - Ai)x, Aiy) = a((A - Ai)x, {Ai - A)y) + a{{A - Ai)x, Ay) 

= a{{A - Ai)x, {Ai - A)y) - b{{A - Ai)x,Ty) 

= a{{A - Ai)x, {Ai - A)y) - b{{A - Ai)x, {T - Ti)y) 

< C5{W,Wi)‘^. 

The second term in (5.4) can be estimated using the following identity 
{St{T - Ti)x, y)H = {Se{T - Ti)x, y - Fiy)^ 
which finally leads to 

\{Si{T — Ti)x, y — Eiy)-u\ < ||5f IjT — Ti\\c(-H)\\I — Pe\\c{H)- 

□ 

The operators T and Ti can be represented by symmetric positive definite 
matrices of dimension N x N (N being the dimension of IT). The following 
theorem is then a standard consequence of matrix perturbation theory (see, for 
instance, [21, Theorem 3, items c) and d)]) and to the equivalences Ai = 1/y-i 
and \i^i = l/yi^i- 

Theorem 5.4. Let J denote the set of indices corresponding to the eigenvalues 
in the cluster W. Then 

sup inf |Ai — Xi j\ < C5{W, Wi)'^. 

i£j iG J 

6 Auxiliary results 

This section contains all technical results which have been used for the proof 
of Theorem 4.1. We arrange the presentation in three subsections: in the first 
one a superconvergence result is proved; in the second one we collect the results 
which hold for all refinements 7/t of a given mesh Fh] finally, in the third one 
we include the results which have been proved for the sequence of meshes {Fi\ 
generated by our adaptive procedure. 
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6.1 A superconvergence resnlt and other useful estimates 

Let II/i denote the orthogonal projection onto Mh- 

Lemma 6.1 (Superconvergence for the source problem). There exist p{h) tend¬ 
ing to zero as h goes to zero such that 

\\Tlnu-TM<pih)\\a-Gf,iT^u)\\^. 


Proof. This result has been proved in [23] and can be found in [29] or [7, §7.4] 
as well. □ 

Let = {1,..., N{h)}\J denote the indices of the discrete eigenvalues not 
belonging to the cluster and assume the initial mesh-size is small enough such 
that 

K = sup sup sup T---r < oo. 

Th jej lA' ~ ^h,k\ 

Lemma 6.2. For all j € we have 

{uh,j,Ti^u) = — - (Tf^u — HhUjUhj)- 

A — Ah,j 


Proof. We have 


-Xh,jiuh,j,T^u) = b{ahj,T^u) = -a{ah,j,Gh(T^u)) = b{Gh{T^u),Uh,j) 
— — X(Tl-hU^ Ufipf 

Adding X{uh,j,T^u) on both sides of this identity leads to 

(A - Xh,j){uh,j,T^u) = X{T^u - T\hU,Uh,j). 


□ 


Lemma 6.3 (cf. [43]). Any eigensolution (A, cr, rt) G K. x S x IT in the cluster 
satisfies 

\\T^u-Ahu\\<K\\IlhU-TM- 

Proof. Let us define Ch = Tjfu—AhU. The expansion in terms of the orthonormal 
basis {uh,j 1 j = 1,..., N{h)} reads 

e/i = ^ otjUh.j with ^ = Ije/ilp. 

j^jc 

This relation, Lemma 6.2, and Parceval’s identity lead to 

^ ^ ^ ^ CtjT 7 {Th'^ ~ ^hUjUhj) 

. _ _ jf-, ^ '^h,j 

/ X 1/2 

□ 
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We are now ready to prove the superconvergence result for the eigenvalue 
problem. 

Proposition 6.4 (Superconvergence for the eigenvalue problem). Any eigen- 
solution (A, cr, u) £ M X E X W in the cluster satisfies 

- A/imII < p{h)\\a - Gh(T^u)\\^. 

Proof. The triangle inequality and Lemma 6.3 give 

WTlhU - Ahu\\ < \\IlhU - T^u\\ + \\T^u - Ahu\\ < (1 + K)\\nhU - T^u\\. 

The right-hand side has been estimated in Lemma 6.1. □ 

The following result contains a useful bound of the norm of the error in E 
in terms of our error quantity. 

Lemma 6.5 (Bound for the E norm). Any eigensolution (A,cr,u) £ R x E x M 
satisfies 


\\a - G^(A?,m)||e < |cr - G,i(A/iu)|a -I- (1 -I- A)||u - A/iu||. (6.1) 

Proof. The stability of the continuous problem implies 

\\(J — G/i(A;iU)||e 

< sup [a{a - Gh{Ahu),T) Pb{(T - Gh{.Ahu),v)+ h{T,u-Ahu)). 
lhlb+lhll=i 

The identity (3.1) together with the continuous and discrete eigenvalue problems 
imply 

6(cr - Gh[Ahu),v) = b{a,v) - b{Gh{Ahu),Ilhv) = X{{P^u,nhv) - {u,v)) 

= KP^u- u,v). 

Estimate (6.1) then follows from the continuity of a and b together with the 
elementary estimate Hm — P^ u|| < ||m — A^mH. 

□ 


6.2 Properties valid for all refinements % of Th 

We start this section by proving the efficiency of our theoretical error estimator 
on the generic mesh 7h. 

Proposition 6.6 (Efficiency). Let {a,u) be an eigenpair associated to the eigen¬ 
value then there exists a positive constant CeS, independent of h, such that 

Ph{u]Th) < Cesd{u,Ahu). ( 6 . 2 ) 

Proof. For the reader’s convenience, we recall the definition of the error indicator 
pLh(u\T) for a given element T € Tn- 

Ph{u]TY = \\hT{Gh{Ahu) - VAhu)\\'^ -b ||/it curlG/,(A,iu)|||. 

+ ^2 ^E\\[GhiAhu)]E ■ tEW's- 

Ee£(T) 
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Following the same arguments as in [10, Lemma 6.3], we can prove that 

h^\\Gh{Ahu) - VAhuW'^ < d{u, Ahuf. (6.4) 

Finally, arguing as in the proof of Theorem 3.1 in [1], we can bound the remain¬ 
ing terms of the error indicator as follows: 

II/iTcurlII^ hE\\[Gh{Ahu)]E-tE\W<\\(J-Gh{Ahu)\\'^^, 

EeSiT) 

where T denotes the union of T and its neighboring elements. 

We then obtain the desired result by summing equations (6.4) and (6.5) over 
each elements T gTh- 

□ 

The next result shows a uniform discrete reliability of the theoretical error 
estimator when evaluated on the mesh Th, refinement of 7h- 

First of all, we recall the well-known discrete Helmholtz decomposition which 
is valid for the hnite element spaces we are considering. Suitable references for 
this result are [2] in the framework of discrete exterior calculus or [33]. In our 
setting the discrete Helmholtz decomposition reads (see [34, Lemma 2.5]): for 
any C.h S S/j there exist ah € Alh and /3h € Vk+iiTh) (the space of continuous 
piecewise polynomial of degree /c -I- 1) such that 

C,h = Ghicth) + curl/3/i. (6.6) 

In particular, ah £ Mh is such that 

a{Gh{cth),Th)+ b{Th,ah) = G \/Th€'S,h , , 

b{Gh{ah),Vh) = b{Ch,Vh) Wvh € Mh- 

By definition of the bilinear form and the fact that div 'Eh = Mh, we have that 
div(G/,(a/j) - (h) = 0, hence Gh{ah) - (h = curl/??;. Using again (3.1) there 
exists Th £ Eh such that divf?; = ah- From the discrete inf-sup condition we 
have ||r/,|] < Clla/il]. Hence 

WcthW'^ = h{Th,ah) = a{Gh{ah),fh) < |G?i(a/i)|a||r?;|| < C'|G?;(a/;)|a||a/t||, 
from which we obtain 

||a/.|| < C|G,;(a/;)U. (6.8) 

Proposition 6.7 (Discrete reliability). Provided the mesh-size of Th is suffi¬ 
ciently small, we have 

\Gh{Ahu) — Gfl-(A/i-it)|a + |[A?iU — Ahu\\ 

< GdreifJ-Hiu-,TH \Th) + G p{H){d{u, Ahu) + d{u,AHu)). 

Proof. From the discrete Helmholtz decomposition (6.6) there exist ah £ Mh 
and Ph £ Vk+iiTh) such that 

Gh{Ahu) - Gh{Ahu) = Gh{ah) + curl/3?;. 
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The term || curl/3?i|| can be bounded by using standard arguments as in [24, 3, 
34]. Actually, taking /3^ as the Scott-Zhang interpolant [40] of Ph on the mesh 
Th, 

I cu.v\Ph\l = a{Gh{Ahu) - GH{AHu),curiPh) 

= -a(G//(A//u),curl(/3/i - /3^)) 

= ^ ( / (/3?i - ;8jj)curlG_y(Ai/M)dx 

TeTH\Th 

-[ {Ph-P h)^h{^hu)- ids). 

JdT ^ 

Standard estimates for the Scott-Zhang interpolant give 

I CVLY\Ph\a < ^lH(u]TH\Th)■ 

The integration by parts and some straightforward algebraic manipulations lead 
to 

\^h{oth)\t = a(G/j(A^u) — Gjy(Aifu), G/i(a/j)) 

= \{P^u-P^u,aH) 

= X{{P^u - n^u, ah) + (n,iU - Uhu, ah - Unah) 

+ {Uhu - P^u,ah)). 

We observe that \\P^u — fl/iit]] < |]A?iU — fljiu]]; indeed, P^u is the best H- 
approximation of u into Wh and is characterized by {P})^u — u, Vh) = {P^u — 
UhU,Vh) = 0 for all Vh £ Wh- Hence, the estimate (6.8), Proposition 6.4, and 
Lemma 6.5 prove for the first and the last term that 

{P^u - UhU, ah) + (Uhu - Ph u, ah) 

< {\\P^u - H^mII + \\Uhu - P^ u]]) lG^{a^)|„ 

< p{H)(d{u, Khu) + d{u, AHu))\Gh{ah)\a- 

For the analysis of the remaining term, set ^ = ah — Unah- It is shown in [34, 
Lemma 2.8 and Equation (3.9)] that ^ satisfies ||^|] < H\Gh{ah)\a- Thus, we 
have with Proposition 6.4 that 

{UhU - Uhu, ^) = {UhU - Ahu, ^) 

= (Hf^w - AhU, ^) + {AhU - Ahu, ^) 

< {p{H)d{u, Ahu) + iL||A;iM - Ai/M||)|Gft(a^,)|a. 

Altogether we obtain for the error in the vector variable that 

|G;i(A;iM) — GHiAHu)\a 

< ph{u; TH\Th) + p{H){d{u, Ahu) + d{u, Ahu)) P H\\AhU - Ahu\\. 

It remains to estimate the error in the scalar variable. 

Let z be the gradient of the solution p of the problem 

Ap = AhU — Ahu in H 
^ = 0 on dU. 
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Using a (non-orthogonal) stable decomposition like the ones adopted in [31, 
Lemma 3.3] or [39, Lemma 2.1], it is possible to find 2 ; G such that 

z = z + curl'i/;. 

In particular we have 

div z = AfiU — Ahu 

ll^:|| + IIV 2 II < \\AhU-AHu\\. 

It follows 

\\AhU - Ahu\\'^ = b{z, AhU - Ahu) 

= b{Ilj^z, Ahu) - billnz, Ahu) 

= -a{Gh{Ahu),Il^z) + a{GHiAHu),Il^z) 

= a{GH{AHu) - Gh{Ahu),Ilf^z) + a(Gi/(A//u), (11^ - 11^)2:) 

< IGniAHu) — G/i(A^M)jQjjnf2:11 

+ a{GH{AHu) — Vh{Ahu), (II^ — n^)^), 

(6.9) 

where we have used the definition of the Fortin operators 11^, IlD, of Gh and 
Gh, and, in the last term, the fact that the quantity a(V//(AijM), (IlD — 11 ^) 2 ;) 
vanishes. 

We observe furthermore that II^z —11^2: = 0 on the unrefined elements Tf/H 
77i. Since 2 : is smooth enough to allow for stability and first-order approximation 
of and nD, we conclude 

jjA/iM — Ahu\\'^ < \Gh{Ahu) — G;i(AhM)jQjjn^2:l| 

+ lliL(G^(Affu) - VH{AHu))\\r^\rjH-Hll^z - Kz)\\ 

< [[A/iU - Ahu\\ 

{\Gh{Ahu) — Gh{Ahu)\a + ^Ih{AhU\ Th \ %))■ 

□ 

By passing to the limit in the statement of Proposition 6.7, and observing 
that for H small enough the second term on the right-hand side can be absorbed, 
we obtain the standard reliability estimate. 

Corollary 6.8 (Reliability). Provided the initial mesh-size is sujficiently fine, 
we have 

^ ^ d{uj, AfiUj) C^ei h-hjUj , Th) ■ 

jeJ j£J 

We conclude this section with a quasi-orthogonality property. 

Proposition 6.9 (Quasi-orthogonality). There exists a constant Cqo such that 
d{AhU, AhuY < d{u, Ahu)"^ — d{u, A/^u)^ -I- Cqo p{h){d{u, AhuY + d{u, Ahu)"^). 
Proof. The proof departs with the following obvious algebraic identities 

\Gh{Ahu) — GHiAHu)\1 = \(T — Gh{Ahu)Y — [cr — Gh{Ahu)\l_ 

— 2a{a — Gh{Ahu), Gh{Ahu) — Gh{Ahu)) 
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\\AhU - AhuW^ = ||u - AhuW^ - ||u - AhuW^ - 2{IlhU - Auu, Auu - Ahu). 

The exact and discrete eigenvalue problems together with the inclusion div Tjh Q 
Mh imply 

a{a - GhiAhu), Gh{Ahu) - Gh{Ahu)) = -b{Gh{Ahu) - (A//u), u - Ahu) 

= X{P^u - Ph u,l\.hu- Ahu). 

Therefore it follows from Proposition 6.4, Lemma 6.5, and the Young in¬ 
equality that 

\a{a - Gh{Ahu),Gh{Ahu) - G_f/(A//u))| -I- |(n?iM - AhU,AhU - Ahu)\ 

< ||n?iU - A/im|| (||A?iM - Ahu\\ + - Phu\\) 

< p{h){d{u, Ahuf + d{u, Ahu)^). 

□ 


6.3 Contraction property 

While the properties of the previous subsection are valid for any refinement 7^ 
of a mesh 7}/, in this section we deal with the mesh sequence 7^ which is the 
output of the adaptive strategy described in Section 3. 

The following property is quite standard and can be proved with the tech¬ 
niques of [17]. 

Lemma 6.10 (Error estimator reduction property). Provided the initial mesh- 
size is sufficiently small sueh that the bulk eriteria for and rp are equivalent 
(see Lemma 3.2), there exist constants pi € (0,1) and K € (0,-|-oo) such that 
T(. and its one-level refinement Te+i generated by AFEM satisfy 

'^pe+i{uj,Ti+if < pi'^pi{uj,Tif + Ki'^d{Ai+iUj, AtUjf'. 
jeJ jeJ jeJ 

The following proposition presents the main contraction property which is 
essential for the convergence proof of the adaptive strategy. 

Proposition 6.11 (Contraction property). Provided the initial mesh-size is 
sufficiently small, there exist p 2 S (0,1) and (3 € (0, -l-oo) such that the term 

= + P^d{uj,AiUjf (6.10) 

j&J 


satisfies 

S.e +1 < P 2 f.£ for all £gN. 

Proof. Throughout the proof, we use the following notation 

e? = '^d{uj,A£Ujf pI = '^pe{uj,ref. 
j&J jeJ 

The error estimator reduction from Lemma 6.10 and the quasi-orthogonality 
from Lemma 6.9 imply the following bound 

P-i+i + 7^1(1 — Cqop{ho))e'j_^i < Pipj + Kiil C'qop(/io))e^. 
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For any e G (0,1), the last bound and the reliability (Corollary 6.8) give 

+ Ki{i C'qop(ho))e^_|_]^ 

< (pi + eC^f,iKi)n£ + iFi(l — e + C'qop(ho))ef. 

We take l3 = Ki{l — Cqop(ho)) and 


P 2 = max 


Pi + eC^gjiFi, 


1 — £ + C'qop(/lo) 
1 Cqo/ 5 (/io) 


SO that 

/if+i + /3e|_|_i < P 2 {p^ + /?ef). 

The choice of a sufficiently small e and of a sufficiently small initial mesh-size 
ho leads to p 2 < 1- □ 

7 Extension to three space dimensions 

The results presented in the previous sections hold true also in three dimensions, 
provided the definitions of the computable and theoretical error indicators are 
modified as follows. 

Definition 9. Let Th be a simplicial decomposition of Ll and let {(Jhp,Uhj) G 
S/i X Mh be a discrete eigensolution computed on the mesh Th- Then, for all 
T G Th we define 

VtijiT)"^ = WhricThj-Vuh,j)\\T + \\hT curl ah,j\\T+ >^F\\[ah,j]F x ufWf, 

FGFiT) 

where hx is the diameter ofT, T{T) denotes the set of faces ofT, hF is the area 
of the face F, and uf is its unit normal vector. As usual, {ah]F x uf denotes 
the jump of the trace of ah x uf for internal faces and the trace for boundary 
faces. 

Definition 10. Let Th G T be a triangulation and let {a,u) be an eigensolution 
associated to the eigenvalue X (in particular, this is used in the definition of 
Ah). For all T G Th we define 

pI{u]T) = |1/it(G(A/iu) - VA,iM)||f. -h ||hTCurlG(A,iM)||f. 

+ ^2 ^F\\[G{Ahu)]F X UfWf- 

FeF(T) 

In the three-dimensional case, the only proof which needs to be modified is 
the one of the discrete reliability of Proposition 6.7 since it relies on the discrete 
Helmholtz decomposition which is different in two or three dimensions. 

Let Vh denote the iL(curl)-conforming edge elements of Nedelec (see [7]). 
Then, in the three dimensional case, the discrete Helmholtz decomposition 
reads (see [34], Lemma 2.6): for any f,h S there exist ah G Mh and fdh G 14 
such that 

(,h = ^h{.oih) + curl ,5/1. 
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Proposition 7.1 (Discrete reliability). Provided the mesh-size of Th is suffi¬ 
ciently small, we have 

|G/i(A?iU) — G/i-(Ajyu)|a + ||A;i(M) — Ah-(m)|| 

< CdreifJ-Hiu;TZ) + Cp{H){d{u, khu) + d{u,AHu)), 

wheren = {T hTh :fnf' 'iT'eiTHXTh)}- 

Proof. Using the discrete Helmholtz decomposition, we write the error in the 
vectorial variable as 

Gh{Ahu) - GniAHu) = Gh{ah) + curl/3?,, 
with ah € Alh and (dh & Vh- 

The term |G?,(Qf?,)|a can be treated without any modification as in the two 
dimensional case. Moreover, following the same argument as in [34, Lemma 3.1.], 
it can be proved that 

I curl fdh\a < pH{u-,n). 

As in the proof of Proposition 6.7, the error in the scalar variable can be 
bounded by using the duality argument of [31, 15] and we can repeat the same 
arguments of the 2D case from Equation (6.9) onwards, concluding the proof. 

□ 

Remark 2. Compared with the two-dimensional case, in the three-dimensional 
version of the discrete reliability, the set Th \Th is replaced with the slighliy 
larger set TZ which essentially is Th \ Th plus one additional layer of simplices. 
The shape-regularity implies that there is a constant C such that 

card(/^) < C car^.{T h \ Th)- 

and therefore the estimate (4.7) remains valid at the expense of the multiplicative 
constant C, and with this modification the proof of Theorem 4-1 applies also to 
the three-dimensional case. 
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